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Abstract. Single pion production in neutrino-nucleon scattering is discussed. The neutrino 
energies characteristic for T2K project are considered. Two new parameterizations of 
axial form factor are proposed. Both of them are obtained in simultaneous fit to ANL and 
BNL data. One of them (which fits better to BNL data) leads to da/dQ^ differential cross 
section which is significantly reduced at low . 



1 Introduction 

Single pion production in neutrino-nucleon/nucleus scattering has been a subject 
of intensive theoretical and experimental studies for many years. One of the first 
comprehensive theoretical model for pion production in electron- and neutrino- nu- 
cleon scattering was proposed by Adler [1]. Since Adler model had appeared many 
other approaches have been proposed. These descriptions were mainly devoted to 
electron-nucleon scattering. Some of them were extended to describe neutrino- 
nucleon interactions. Let us mention two of the newest approaches: Sato and Lee 
model [2] and the model proposed by Hernandez et al. [3]. Both of them are applied 
to describe charged and neutral current neutrino-nucleon interactions. 

From the theoretical point of view there are useful analogies between electron- 
nucleon and neutrino-nucleon interactions. The vector part of charged electro-weak 
hadronic current is related with the electromagnetic hadronic current - it is a result 
of conserved vector current (CVC) hypothesis. For electron-nucleon scattering high 
accuracy experimental data are available. This allows to verify and fine tune the- 
oretical approaches. Well tested models in electron-nucleon interaction should de- 
scribe well the vector contribution to i/-nucleon scattering, but for neutrino-nucleon 
scattering more important is the proper modeling of an axial contribution. The de- 
scription of the axial contribution to the cross section can be mostly verified by 
comparing with neutrino-nucleon scattering data. 

Neutrino-matter interactions have been intensively studied for many years but 
experimental data for Itt production are limited. Most of the data sets have been 
collected for charged current vN scattering, and there are only a few total cross 
section points for neutral cuiTent neutrino-nucleon scattering [4] [5]. A study of the 
axial contribution, a fitting of the axial form factors, requires analysis of da/dQ^ 
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differential cross sections. Such data have been collected by several experiments. 
Let us mention two of the most important experiments, namely the bubble chamber 
experiments, ANL [6] (Argone National Laboratory) and BNL (Brookhaven Na- 
tional Laboratory) [7]. 

In the ANL experiment the differential cross section da/dQ^ was measured up 
to = 1 GeV^ while in BNL experiment data were collected up to = 3 GeV^. 
In both projects deuterium was used as a target and very similar measurement tech- 
niques were applied, but the collected cross section data seem to be in a disagree- 
ment. The BNL total cross sections are systematically higher (by 20%) than the 
ANL total cross sections [8]. Additionally fitting to ANL da/dQ^ data leads to a 
functional form of the axial form factor which has different dependence than the 
one extracted from the BNL data [9]. Therefore, usually (in the literature) the axial 
form factor fits to either the ANL or BNL data. 

During the last few years the experiments K2K (KEK to Super-Kamiokande 
long -baseline neutrino oscillation experiment) [10] and MiniBooNE [11] (Booster 
Neutrino Experiment at Fermilab) have collected new data also for single pion pro- 
duction. Since in K2K a water cherenkov detector (Super-Kamiokande) was used 
neutrino-oxygen interactions were mainly detected, while at MiniBooNE neutrinos 
mainly interact with carbon nuclei. Modeling of neutrino-nucleus scattering requires 
to consider some description of the nucleus, which must be implemented in the data 
analysis. 

K2K and MiniBooNE experiments are different in the design and they used dif- 
ferent measurement techniques but in both projects it has been observed that the 
da/dQ^ data points at low (four-momentum transfer) are below theoretical pre- 
dictions (Monte Carlo (MC) simulation) [12] [13]. This discrepancy between theory 
and experiment can be explained by assuming that the contribution from the coher- 
ent pioi{3 production is negligible [12]. On the other hand it is possible that this 
effect can be explained by more careful discussion of nuclear effects. There are 
also possible other sources of this disagreement like improper account of the lep- 
ton mass effects [14]. Eventually, the dependence of the cross section is mainly 
constrained by the axial form factor 

The investigation of Itt production in neutrino-nucleon scattering is important 
to understand the character of neutrino-matter interactions as well as the structure 
of the nucleons and nuclei. The subject is also important from the practical point of 
view: a proper prediction of cross sections for Itt production plays a crucial role in 
data analysis of long baseline oscillation experiments. Indeed in K2K the oscillation 
of neutrinos was investigated by observing distortion of the energy spectrum in 
the far detector [15]. The energy spectrum was reconstructed by considering quasi- 
elastic (QE) charged current events but some processes with single pion production 
stand background to these events. 

In few years T2K [16] (Tokai to Kamioka), a next generation long baseline neu- 
trino oscillation experiment, will start collecting data. A study of i/^ ~^ oscil- 

^ In the coherent interaction the quantum numbers of the nucleon are unchanged and pions 
are produced mostly in the forward directions. 
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lation is the main goal of the project. It can be only investigated by observing the 
appearance of the electron neutrinos in the far detector In practice, to observe elec- 
tron neutrinos the electrons produced in charged current j^e-target interactions must 
be observed. Simultaneously the detector will be able to detect vr^s which are pro- 
duced in neutral current neutrino-target interactions. This is important background 
for i/g measurement because 7r° decays into two photons which can be misled with 
the electron shower Since a small number of electron neutrinos is expected in the far 
detector, success of the experiment may depend on proper estimation of the number 
of 7r°s produced in neutral current neutrino interactions. 

The description of tt*' production in neutral current neutrino-nucleon scattering 
still requires careful attention. As has been akeady mentioned there are limited ex- 
perimental data for these processes. Standard model gives a hint on how to construct 
the hadronic current for neutral current interaction. The neutral current (NC) is ex- 
pressed thought the third isovector component of the vector-axial current and the 
electromagetic current 

J^C ^jCC,U_^^^^2q^jEM^ (1) 

where Qy^ is the Weinberg angle. The relative weight of the vector and axial con- 
tributions to the neutral current cross sections is different than in the case of the 
charged current reactions. One can imagine a model which fits well to charged cur- 
rent data, but overestimates axial and underestimates vector contributions. In that 
case predictions for neutral current reactions will be wrong. 

In this talk we will present comparison of two descriptions for single pion pro- 
duction in neutrino-nucleon scattering, namely the Rein Sehgal (RS) model [17] and 
the isobar formalism [18]. The first model is implemented in the Monte Carlo codes 
which were/are used in data analysis of K2K and MiniBooNE experiments. We will 
discuss how to improve the Rein Sehgal model to get more precise description of the 
charged current and neutral current cross sections (in particular for vr" production). 
Our idea is to consider new vector and axial form factors. We express the RS model 
form factors by the ones from the isobar formalism. The new vector contribution 
fits well to the electoroproduction data, while the axial contribution was obtained in 
simultaneous fit to the ANL and BNL differential cross section da/dQ^ data. 

2 Single pion production 

One can distinguish two mechanisms for single pion production in neutrino-nucleon 
scattering: 

(i) resonant - the nucleus is excited to the resonance state, then the resonance de- 
cays to a single pion and nucleon. It is illustrated in Fig.[T^); 

(ii) nonresonant - the pion is produced without resonant excitation of the nucleon. 
It can be described by Born diagrams (Fig.lTJ))) [19]. 

Since the resonance is not directly observed (only its decay products), it is impossi- 
ble to distinguish between both mechanisms. 
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Figure 1. In Fig. (a) Feynman diagram for the resonance pion production in uN scattering is 
sliown. In Fig. (b) one of tiie typical Born diagrams for the nonresonant pion production in 
vN scattering is presented. 

However, it seems that for the reaction: 



the nonresonant contribution is neghgible - to fit the hadronic invariant mass dis- 
tribution of the TT+p system it is enough to consider the Breit-Wigner distribution 
of the P33 (1232 resonance [7]. Therefore process (|2]l seems to be appropriate for 
fine tuning the description of the resonance contribution to the pion production in 
i^N scattering. For other possible channels the nonresonant contribution as well as 
heavier than 1232 MeV resonances must be taken into account. However, if the neu- 
trino has energy around 0.7 GeV (characteristic energy for T2K project) it will be 
enough to consider only the P33(1232) resonance. 

In this talk we present results of comparison between the RS model and the iso- 
bar formalism for the P33 (1232) excitation. Both descriptions should lead to similar 
resonance helicity amplitudes for the P33(1232) excitation. 

2.1 Rein Sehgal model 

The Rein Sehgal model is an extension of the relativistic harmonic oscillator quark 
model (FKR - Feynman, Kislinger and Ravndal) [20] . The FKR model was applied 
to describe photo- and electro- production of resonances [21], and then neutrino- 
production of resonances [22], [23]. 

In the FKR model baryon wave functions are given by symmetric representa- 
tions of SUfiavorC^) X SUspin{2) X 0(3). The nucleon and Z\(1232) are ground 
states of the harmonic oscillator. Other excited states of the nucleon are classified 
by: level in the harmonic oscillator model, spin, orbital angular momentum etc.. 
Feynman et al. proposed an operational form of one body quark currents for elec- 
tromagnetic and weak transition of nucleon to resonance state. Knowing the nucleon 
and resonance wave functions one can compute the following helicity amplitudes: 

It is the 3/2 spin resonance denoted also by Zi(1232). 



v + p 



^J. +Tr^ +p. 



(2) 
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h = {N*,s'\Jt + —J,\N,s), (3) 

Qres 

f± = {N\s'±l\J±\N,s), J± = T^iJx±iJv), (4) 

where = — g^, = (i/res, 0, 0, g^es) and Ures, qres are energy, momentum 
transfers computed in the resonance rest frame, s and s' denote spins of incoming 
and outgoing particles. 

The hadronic charged electro- weak current has vector - axial structure: J^^ = 
— J^. In the formalism proposed by Feyrmian et al. the vector and axial one 
body quark currents are multipUed by two phenomenological functions: the vector 
and axial form factors 

J — > (^v,aJ 

Final hadronic state can be either a nucleon or a resonance. It means that formalism 
should be able to reconstruct helicity amplitudes for elastic ep scattering and quasi- 
elastic m scattering. This property was applied to find expressions for Gv,a- They 
are proportional to the elastic vector and axial nucleon form factors: 

o.w^) = (i + :|)"(i + 3S.)*. 

where My and Ma are vector and axial masses while M denotes nucleon mass. 

Rein and Sehgal extended the FKR model to describe single pion production 
induced by charged and neutral currents neutrino-nucleon interactions. Pions are 
produced by resonance excitations but the description of the nonresonant contribu- 
tion was also proposed. In the RS approach 18 resonances with hadronic invariant 
mass smaller than 2 GeV are discussed. 

It must be emphasized that the RS model successfully described neutrino cross 
section data if some reasonable accuracy was assumed. However, higher precision 
measurements require more detailed description of the cross sections. We noticed 
that the RS model underestimates inclusive F2 ep structure function [9] and con- 
sequently ep cross sections. Therefore we proposed an effective description of the 
vector contribution, which more accurately fits to the ep experimental data. 

As was mentioned above next generation long baseUne experiments wiU mea- 
sure neutrino-nucleus scattering mainly in the Z\(1232) resonance region. Therefore 
our improvements concern the ^33(1232) resonance. 

2.2 Isobar formalism 

In the FKR model the axial and vector form factors were obtained by referring 
formalism to elastic and quasi-elastic scattering. But one can propose a different idea 
for introducing form factors into the model. Since the zi(1232) resonance region is 
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(6) 
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Figure 2. da / dQ^ differential cross sections for reaction + p — > + 7r+ + p. In the left 
figure cross sections computed for the ANL beam are presented [6]. In the right figure the 
differential cross sections are computed for i5 — 0.7 GeV - an average neutrino energy for 
T2K experiment. The cross sections are computed with the original RS form factors (dotted 
lines), and our fits. The solid and dotted lines denote cross sections obtained with l ll7t and 
dl6t axial form factors respectively. The cut on the invariant hadronic mass W < lA GeV is 
imposed. 



of our interest, it is natural to use knowledge about the P33(1232) production as 
reference point for the FKR/RS form factors. For that reason it is useful to compare 
the isobar formalism for the P33(1232) excitation [18] with the RS approach. 

The isobar formalism gives very phenomenological, only in terms of form fac- 
tors, description for the Z\(1232) excitation. Helicity amplitudes computed with this 
approach can be easily compared with corresponding once in the RS model. Since 
the isobar model has been updated [24] lastly, the comparison may show a hint on 
how to improve the RS model predictions. 

In the isobar formalism the hadronic current for the P33(1232) excitation is 
constructed to satisfy Lorentz invariance [18]. The initial nucleon state is given by 
Dirac spinor while the P33(1232) state is the Rarita-Schwinger field. The vector 
current is conserved and satisfies CVC hypothesis while axial current is constrained 
by PCAC hypothesis. These assumptions reduce the number of independent form 
factors. 

The vector current for (1232) excitation is expressed by three form factors 
CY, CY and 

{A++ip')\ \N{p)) = ViMp') Kt^q" " q^T^] 75uip), (7) 

where 

fjV (jV 

= ^7. + + ^P., (8) 

and g^^ is Minkowski metric. 



Resonance contribution to single pion production in neutrino-nucleon scattering 



7 




E [GeV] E [GeV] 

Figure 3. Total cross sections for two neutral current reactions: u + n u + n" + n (left 
figure) and + n ^ u + n~ + p (right figure). The solid and dashed lines correspond to 
cross sections computed with l ll7t and M6\ axial form factors respectively, while by dotted 
lines cross sections obtained with the RS form factors {Ma = 1.1 GeV) are shown. The data 
is taken from Ref. [4] (black squares) and Ref. [5] (white squares). The cut on the invariant 
hadronic mass W < 2 GeV is imposed. 



Axial current depends also on three form factors C4 , C5 and Cg 



{A++ip')\j^f\Nip))=V3^xip') 
where 



From the PCAC hypothesis 



uip), 
(9) 

(10) 



where m^r is pion mass. Usually the Adler relation = —C^ /A is also assumed. 
Then the axial current depends only on C^{Q^) form factor. 



3 Comparison of RS model with isobar formalism 

To compare the RS model with isobar formalism three independent helicity ampli- 
tudes (see expressions: |3]|4]l are evaluated in both approaches. It is well known that 
the quark model is not able to reproduce a charged contribution of the ^33(1232) 
resonance (expression ^ vanishes). Therefore, the exact comparison is impossible 
but it is still reasonable because the contribution from helicity amplitude (O is very 
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small. Comparing helicity amplitudes allows us to express Gy by , CX and 
functions (for more details see Ref. [9]): 



Gl'^{W,Q^) = y3(Gi%W,Q^)y + {G(^iW,Q^))\ (11) 



where 



2^3 



2M2 



2M2 



Gl\W,Q^) = - 



F{Q^W) 

'3 



2^3 



2M2 



= — O4 
^ 2M2 



+ Q2 _ ^2 _^ _^ g2 
G3 - 



MW 



F{Q\W)=[1 + 



(M + Wf 



(12) 

• (13) 
(14) 



where by W hadronic invariant mass is denoted. 

Lalakulich et al. [24] did some effort to fit , CX and CX form factors to 
resonance photo-production data. The obtained parameterizations lead to a vector 
contribution which agrees with elastic ep experimental data, and with MAID predic- 
tions [25]. It seems that a use of these form factors guarantees adequate description 
of the resonance vector contribution in neutrino-nucleon scattering. Therefore we 
express Gy form factor by Lalakulich et al. parameterizations. 

Similarly as in the case of vector contribution three independent hehcity am- 
plitudes are computed for axial current. Here, the relations between amplitudes are 
more complicated than in the vector current case and it was possible only to find an 
approximate relation between Ga and C^: 



\W,Q^) = ^F{Q\W) 



1- 



8M2 



(15) 



where, in our notation Ga = (0.76) • Ga- 

To get a proper form of the (Q^) axial form factor we proposed a numerical 
procedure [9] which allows us to find simultaneous fit to both ANL and BNL data 
sets. We fixed ANL cross sections normahzation then the normahzation of the BNL 
data had to be changed in order to obtain a consistent fit. 

We discussed two different parameterizations of the axial form factor. 

In the first case we assumed that (0) is constrained by PCAC, and equals 1 .2. 
For that case we obtained: 
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with axial mass ==0.54 GeV^. 



We consider also parametrization with C5 (0) treated as a fit parameter. In that 



case we assumed that has functional form: 



c^iQ')-- — (17) 



We got C^(0) w 0.88, with Af ^ w 9.71 GeV^ and w 0.35 GeV^. It is interest- 
ing to notice that the value of at = is very similar to (0) = 0.867 which 
was computed in more advanced approach (chiral constituent quark model) [26]. 



4 Numerical results 

In this section we show how the new form factors, described in the previous section, 
change differential and total cross sections for neutrino-nucleon scattering. 

In Fig. I2] we plot cross sections computed with new form factors and the RS 
model form factors. In the left panel of Fig. |2] theoretical predictions are compared 
with da/dQ^ ANL data for reaction v + p ^ /i~ + 7r+ + p. The cross section 
computed with the axial form factor ( fT6b fits better to ANL data than the fit given 
by Eq. Small value of axial form factor O at = q (C^(0) = 0.88) leads 
to significant reduction of da/dQ^ cross section at low Q^. Let us also remark that 
this form factor fits better to the BNL da/dQ'^ data (see Fig. 3 in Ref. [9]). 

In the right panel of Fig. |2] we show the differential cross sections computed 
for neutrino energies E=0.7 GeV - an expected averaged energy of T2K neutrino 
beam. Similarly as in the left panel, we present cross sections for charged current 
reaction 1/ + p ^ /i^ + 7r+ + p. We see that three different parameterizations of 
the axial form factor lead to different dependence of differential cross section. 
Let us notice that the the original the RS axial form factor is a kind of compromise 
between the axial form factors preferred by ANL and BNL data. From that point 
of view it seems that a use of the Rein Sehgal description leads to world averaged 
values of neutrino cross sections. 

It has been already mentioned that the cross sections for pion production in 
neutral current neutrino-nucleon scattering are of extreme interest. In Fig. [3] we plot 
total cross sections for two neutral current reactions: v + n r/ + 7r" + n and 
V + n ^ ly + + p. The cross sections computed with the form factors presented 
in this talk are by 20% smaller than predictions obtained with the original RS model 
(with axial mass LI GeV). However, the small number of data points did not allow 
to verify which approach fits better to the data. 



5 Final remarks 



In this talk we presented some improvements of the RS model. We proposed modi- 
fications of the RS description which correct description of the vector and axial con- 
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tributions. Our improvements can be easily applied to existing Monte Carlo codes. 
We believe that after our modifications the RS model more precisely describes total 
and differential cross sections. 

We showed also that it is possible to find a simultaneous fit to both ANL and 
BNL data, but the normalization of BNL data had to be changed by about 25%. 

One of our parameterizations of the axial form factor iT7\ leads to the da/dQ^ 
differential cross section which is significantly reduced at small Q^. It is the 
parametrization preferred by BNL data. It is possible that a use of this form factor 
in the data analysis may help to understand the observed at K2K and MiniBooNE 
discrepancy between experimental data and theoretical predictions at low Q^. 
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